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ABSTRACT 



Propagation of extremely short unipolar pulses of electromagnetic field 
( "videopulses" ) is considered in the framework of a model in which the mate- 
rial medium is represented by anharmonic oscillators (approximating bound 
electrons) with quadratic and cubic nonlinearities. Two families of exact ana- 
lytical solutions (with positive or negative polarity) are found for the moving 
solitary pulses. Direct simulations demonstrate that the pulses are very ro- 
bust against perturbations. Two unipolar pulses collide nearly elastically, 
while collisions between pulses with opposite polarities and a small relative 
velocity are inelastic, leading to emission of radiation and generation of a 
small-amplitude additional pulse. 

PACS number: 42.65.Tg 
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1 INTRODUCTION 



Very short nonlinear pulses of the electromagnetic field, which contain few 
optical cycles |], ||, p|, £|, || , or even of half a cycle [f| , have recently attracted 
a great deal of attention. They are referred to as ultrashort pulses (USP), 
or videopulses |l|, |j, |J. From the viewpoint of applications, the interest to 
these pulses is stipulated by a perspective of a drastic increase in the rates 
of data transmission and procession. One of possibilities for generation of 
USPs is compression of a preliminarily phase-modulated pulse. The phase 
modulation can be induced by the Kerr (cubic) nonlinearity . For example, 
pulses with the duration r p = 4.5 fs, consisting of two optical cycles with the 
frequencies belonging to the visible range, and with r p = 40 fs, containing 
just one optical cycle at the frequency of the CO2 laser, were reported in Ref. 

!■ 

If the width of a pulse is much larger than the period of the optical 
cycle, the pulse can be naturally described by means of its slowly varying 
envelope, whose evolution is governed by an effective parabolic (nonlinear 
Schrodinger) equation J7J. However, this approximation cannot be applied 
to USPs, which calls for development of other theoretical models admitting 
a relatively simple description of the propagation of extremely short pulses. 

The existing models that can be used for the description of USPs fall into 
two classes, resonant and nonresonant ones. A resonant medium is modeled 
by an ensemble of atoms with discrete energy levels. The field frequency 
being close to a frequency of an atomic transition, only the (nearly) resonant 
component of the field is taken into account. In particular, the simplest 
model of the self-induced transparency assumes a medium consisting of two- 
level atoms resonantly coupled to the field. In the case of multi-level atoms, 
it may be necessary to consider several monochromatic components of the 
field with different frequencies. 

In the works || [11]] , USP solutions to a closed system of the Maxwell- 
Bloch equations have been found analytically without assuming a separation 
into the carrier wave and envelope. Generally, the Maxwell's equations ad- 
mit the propagation of electromagnetic waves in both directions. If, however, 
the nonlinear contribution into the medium's polarization is small, the uni- 
directional wave propagation may be assumed. This approximation admits 
reducing the wave equation to the first-order one without any assumption 
about the shape of the waves. 

The nonlinear dynamics of the medium driven by the electromagnetic 
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filed is frequently modeled in terms of anharmonic oscillators. In particular, 
the propagation of a linearly polarized USP was considered in Refs. 0, [12| 
in the framework of the material model based on the Duffing oscillators, 
the nonlinear response of the medium being cubic. The propagation of a 
linearly polarized USP in a dispersive medium was modeled by means of 
quadratically nonlinear oscillators in Ref. |l| (without dispersion, this was 
done in Ref.0). 

In all these works, the oscillator represents a response of the high-frequency 
electron degree of freedom to the electromagnetic field. Since the USP spec- 
trum is concentrated at low frequencies, ion oscillations may also give a con- 
siderable contribution to the full polarization of the medium. Nevertheless, 
we will neglect the ion component in the material response (the propagation 
of femtosecond pulses in a medium with the nonlinearity determined by both 
electronic and ionic (Raman-scattering) degrees of freedom was considered 
in Ref. 0). 

An objective of the present work is to study the unidirectional propa- 
gation and interactions of linearly polarized USPs in a nonlinear dispersive 
medium modeled by an anharmonic oscillator combining quadratic and cubic 
nonlinearities. As is well known, in the case when the oscillator is quasi- 
harmonic, the quadratic and cubic nonlinear terms produce effects of the 
same order of magnitude [JT5], that is why it is natural to consider a mixed 
model of this type. It is also relevant to mention that modeling dynamics of 
the broad (rather than ultrashort) optical solitons in a medium with com- 
peting quadratic and cubic nonlinearities has recently attracted considerable 
attention, see, e.g., the works Jl6| and references therein. 

The rest of the paper is structured as follows. The model is derived in 
section 2. Two families of moving USP solutions with positive and negative 
polarities are found analytically in section 3 (only one family survives in the 
limit of a vanishing cubic nonlinearity). Stability of the pulses and their 
collisions (for both signs of the relative polarity of the colliding pulses) are 
investigated in section 4 by means of direct simulations. 



2 THE MODEL 

The one-dimensional propagation of the electromagnetic waves in a nonlinear 
medium is governed by the wave equation, 
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d 2 E 1 8 2 E _ 4vr d 2 P 

Ih 2 ~ ~ ^9^' (1) 

where P is the polarization of the medium. In the framework of the usual 
unidirectional approximation [TJ|, Eq. (|I|) can be reduced to the first- 
order equation, 

dE 1 dE _ 2tt dP 

~dz~ + c~dt~~~~dt' ^ 
We adopt a simple anharmonic-oscillator model for the medium, which 
is a frequently used approximation fl9j (see also |20], If X as the dis- 



placement of an electron from its equilibrium position, the motion equation 
(which neglects friction) can be written as 

^ + u 2 X - k 2 X 2 + k 3 X 3 = 2-E(z, t), (3) 

where u Q is an eigenfrequency of the oscillator, while k 2 and k% are anhar- 
monicity coefficients. The term on the right-hand side of Eq. repre- 
sents the force exerted on the electron by the electromagnetic field, where 
£ = (e + 2)/3 is the Lorentz factor, e is the electron's electric charge, and 
m is its mass. However, it is possible to absorb £ into an effective mass 
m e g = m/£ . Hereafter, we will use m as a symbol for this effective mass. 
Finally, the dynamical variable X is related to the medium's polarization, 
P = n^eX, where is the density of the oscillators (atoms). 
We define rescaled variables, 

C = z/l, r = iu (t- z/c), £ = E/A , and q = X/X , (4) 

where A = m^o/el^l, Xq = Uq/\k2\, 1// = 2nnAe 2 / (mcuo) = u 2 /2u c, 
and 

u p = (A7in A e 2 /m)- 1/2 (5) 

is the plasma frequency. In terms of the rescaled variables, Eqs. (0) and (|3]) 
take the form 

d£ dq d 2 q 2 3 

dC=~^ d^ 2+q - q +2M =s > (6) 

with the single remaining parameter /i = (kzuJq/2k%)- Equations (§) furnish 
a final form of the model. 
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To estimate the magnitude of parameter /x, we consider Morse potential 



U M (X) = U [exp(-2X/r ) - 2 exp(-X/r )], 



where ro is an atomic characteristic size and Uq is a ionization potential. 
The expansion of this potential in terms of a power series in X leads to 



= 2U /mr%, \k 2 \ = 3U /mr%, k 3 = 7U /3mrl Thus \i = 7/27 « 0.259. 
Consider now Kepler's potential with Coulomb long-range attraction: 



U K {X) = u 



r 



r + X 



r 



r + X 



In this case we have Uq 



where r corresponds to equilibrium position. 
2C/o/wtq, \kq\ = 6Uo/mrQ,K 3 = YlUo/mrQ , and \i = 1/3. In following 
numerical simulation of the steady state pulses propagation about the same 
magnitude of parameter \i will be used. 



3 DYNAMICAL INVARIANTS 

It is worthy to note that the system of equations @ can be derived as the 
Euler-Lagrange equations from the action functional 

S — J C[q, 4>}dxdt, 

where the Lagrangian density is 

1<90 <90 1 (dq\ 2 1 2 1 3 n 4 d<p 
C= 2^m + 2{^) ~2 q + 2 q -~2 q +q Tx (7) 

and the independent variables were temporarily (up to the end of this section) 
redenoted as ( — > t, r — > x, in order to have the notation similar to that in 
the classical field theory (a definition of </> is given below). 

Application of the variational procedure to the action S yields equations 

d 2 <p d 2 q d 2 q 2 3 d(f) 



dtdx dx 2 ' dt 2 11 ^ q dx 
Identifying as a potential for the fields q and S, so that q = —d<j)/dt and 
£ = d(j)/dx, makes these equations identical to Eqs. (J6|), which can be further 
transformed into a single equation 
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Note that Eq. (|9]) may be regarded as a continuity equation, 

!(,_, + W + g) + |_a do, 

Integration of Eq. ( |T0D in rc leads to 

+00 

|/( 5 -^ + 2M 3 )dl + |(| K ) +5K = . 

— OO 

Assuming zero boundary conditions for q and its derivatives at x — > ±00, 
the latter relation implies the conservation of a dynamical invariant, 

+00 

h= J (q-q 2 + 2fiq 3 )dx. (11) 

—00 

By taking into account the second equation from Eqs. (|^) in the form 

+ g-g 2 + 2^ 3 = £, (12) 
one can find that l\ is the "area of the pulse" , 

/+00 
Sdx. (13) 
-00 

To find another dynamical invariant, we multiply Eq. (||) by q and trans- 
form the resulting relation into a form 



dt 



1 2 2 3 3/4 4 I (dq 
2 q ~3 q + Y q -2U 



d_ A 2 <9 2 g \ 
cte \2^ q dxdt J 



0. 



(14) 



Thus, a new continuity equation is obtained, giving rise to a new dynamical 
invariant, 
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1 9 2 o 3u , 1 / dq 



2\dx 



dx. 



(15) 



Let us now get back to the conserved density in the continuity equation 



(H), 



Q = g-g 2 + 2/ig 3 + 0. 



Multiplying both sides of Eq. (10) by Q, we obtain 



d (\ 



dt \2 



Q z )+Q 



dq 

dx 



This immediately leads to a continuity equation 



d (\ 



-,Q 2 + 



d_ 

dx 



1 3 , V- 4 . 1 ( d( l 



dt \2 

and the third dynamical invariant 

+oo 



1 



2 



(16) 



(17) 



Taking into account the relation (|T2"|), this integral may be interpreted as a 
"pulse energy", 



(18) 



£ 2 dx. 



To find the interpretation of the dynamical invariant given by Eq. (]T5j) , 
we resort to the Lagrangian density (^). The density of the canonical Hamil- 
tonian for this dynamical system can be obtained from £ by means of the 
standard Legendre transformation, 



U 



dC 



dC 

dq,t 



q,t - £ 



1 (dq 

2 1 dx 



1 2 1 

H — q q 

2 H 3 y 2 



2 x 3 , A 4 4 c 

9 ~~ OP- 



The variable £ can be eliminated from it, using Eq. (|12"D , so that 
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"--sC-SH®'-^-^' (19) 

Omitting the full derivative, the Hamiltonian corresponding to the density 
( fl9| ) takes the form 



H 



1 /9g 

2 V9x 



dx = — In 



cf. Eq. flT5|). Thus, the dynamical invariant — 12 is nothing else but the 
canonical Hamiltonian of the system under consideration. 



4 ANALYTICAL SOLUTIONS FOR THE UL- 
TRASHORT PULSES 

It seems plausible that the system of Eqs. @ is not an integrable one. 
Nevertheless, some exact analytical solutions, describing the propagation of 
USPs, can be found. To this end, one assumes that £ and q depend on a 
single variable, 

V = T-<Z/a = u (t-z/V), (20) 

with some constant a. An expression for the velocity V of a steadily moving 
pulse then follows from Eq. (§]), 

1 1 

V ~ c 

Up being the plasma frequency defined by Eq. (H) . The first equation of the 
system (|6]) can be integrated to yield 

£ = aq. (22) 
Next, the second equation from the system (^) takes the form 

_ ( Q _ i) q _ g 2 + 2fiq 3 = 0> (23) 
This equation can be integrated once, 



2a \uj J 
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— I -(a-l)g - -g + /xg = const, 



(24) 



As we are interested in solitary-wave solutions, it is necessary to set const 
0, so that (E|) yields 



(25) 



A substitution q = 1/y transforms Eq. (p5|) into 



dr] 



\ 



(a-1) 



1 



V 



l + 9//(a- 1) 



3(a-l)7 9(a-l) 2 
A new dependent variable £ will be used, denned by 

1 



(26) 



y + 



3{a- 1, 
with a — ±1. In terms of £, Eq 



9(a- l) 2 
is 



1/2 



e„ 1 + 9Ma ~ 1) ) 



from where one finds £(77) = cos1i(a/ (a — 1)77) . 

Thus, we have obtained a family of exact solutions parameterized by the 
continuous positive parameter (a — 1) and discrete one a = ±1, 



g (+) (r?; a) 



3(a- r 



a/1 + 9(a - l)/xcosh(^/(a - 1)77) - 1 ' 
3(a- 1) 



(27) 
(28) 



'1 + 9(a - l)/xcosh(y / (a - 1)77) + 1 ' 

the superscript standing for a. 

The pulses represented by the solutions ([27]) and fl2"%D have different polar- 
ities. In the limit /x — > 0, corresponding to the model with a purely quadratic 
nonlinearity, Eq. (|27|) goes over into a singular solution, 

3(a - 1) 



2sinh 2 (v / (a- 1)77/2) 
10 



while Eq. (|28D yields a nonsingular pulse in the same limit, which was already 
found in Ref. f23fl , 



(-) 3(a 



9quadx(«^) 



2cosh 2 (^(a- 1)77/2) 



Thus, in compliance with the results of Ref. |23| , only one family of nonsin- 
gular pulses exists in the case of the purely quadratic nonlinearity. 

By using Eq. (^3), we retrieve an expression for the electromagnetic field 
corresponding to the solutions (p7|) and fl28|), 



gW( a - n) = 3a(a-l) 

v/l + 9(a - l)/icosh(v/(a - 1)77) - 1 ' 

5H(a; ,) = - 3a{a ~ l \ . (30) 

y/l + 9(a - ^//cosh^a - 1)77) + 1 

We stress that the velocities of two pulses which have the opposite polar- 
ities but the same value of (a — 1) are equal, since they are determined by 
Eq. (|SID and depend only on a. 



5 NUMERICAL SIMULATION OF PROP- 
AGATION AND COLLISIONS OF THE 
PULSES 

Although the pulse solutions have been obtained in the exact form, their sta- 
bility and interactions should be studied by means of numerical simulations. 
To simulate the system (|5]) , it is convenient to transform it into the following 
form: 

I = -* < 31 > 
| = p , ! = f -, + ,*- 2/ y. 02) 

Since Eq. ( PT| ) contains only the derivative with respect to (, and the deriva- 



tives with respect to r are present only in Eqs. (|3~2"D , one can use any tech- 



nique of numerical integration of ordinary differential equations. For solving 
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Eqs. (P2"|), we employed the prediction-correction method. The fourth-order 
Runge-Kutta routine was used for integrating Eq. ([H]). As the initial con- 
ditions, we took the analytical solutions given by Eqs. ( p7|) and (|29]), or ( |28| ) 
and (|30p, at ( = (with regard to Eq. (|20|)). The boundary conditions are 
q(C, t) = P(C, t) = at \t\ — ► oo. 

By using the Morse Um(X) and Kepler Uk{X) potentials one can es- 
timate the normalized amplitudes A and X . For Morse potential A = 
(4/3)(f/o/er ), X = (2/3)r , and for Kepler's potential we have A = (2/3)(f/o/er ) 
(l/3)r . In both case (Uo/er ) can be identified with atomic field pl| . Con- 
sequently, the peak amplitude of the steady state pulse must be not exceed 
(Uo/er Q ). Furthermore, the anharmonic-oscillator model proposes that the 
magnitude of q is less than one. It leads to the inequality (a — 1) < 1. How- 
ever, in order to obtain the descriptive illustrations of the numerical results 
the parameter a have been taken both from the interval 1.01 < a < 1.1 and 
from the interval 2 < a < 8. 

A typical numerical result, with a = 7.25 and /i = 2/9, is displayed in 
Fig.l, which demonstrates stable propagation of the pulse, at least up to 
C = 200. To further test the robustness of USPs, we added to the initial 
configuration perturbations in the form of a quasi-harmonic or biharmonic 
wave with amplitudes < 0.1 of the pulse's amplitude, and with different fre- 
quencies. As a typical example, Fig. 2 demonstrates evolution of the pulse 
with a = 5 and fi = 1/3 , with the initially added perturbation in the form of 
a long quasi- monochromatic packet filled by the wave SS(t) = <5£o sin(15r). 
As is seen in from Fig. 2, USP is not scathed by the perturbation. With the 
increase of the amplitude of the quasi-harmonic perturbation, we observed no 
conspicuous change of USP. Simulations of the pulses evolution in the pres- 
ence of a biharmonic (two-frequency) perturbation (not shown here) showed 
that the increase of the amplitudes of the biharmonic perturbation results in 
a small change of the value of the parameter a corresponding to the finally 
established USP. A similar effect was generated by an initial perturbation 
in the form of an additional rectangular pulse superimposed on USP, see an 
example in Fig. 3. The increase of the perturbation's amplitude leads to an 
increase of the final value of the USP's amplitude and decrease of its velocity. 

Simulations of interactions between pulses with widely different values of 
the parameter a (i.e., velocities) have shown that their collisions are quasi- 
elastic, irrespective of the polarities of the colliding pulses: after passing 
through each other, the pulses retrieve the same shapes and velocities as they 
had before the collision, the only result being a shift of their centers. How- 
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ever, the character of the interaction between pulses with equal and opposite 
polarities becomes different with the decrease of their relative velocity. In the 
case of equal polarities, the pulses with a small relative velocity demonstrate 
strong mutual repulsion: the distance between them attains a minimum, 
and then they start to separate again, so that they never completely over- 
lap. Strong energy exchange between the pulses takes place around the point 
where they attain the minimum separation. The energy exchange gives rise 
to mutual interconversion of the two pulses, so that after the collision they, 
effectively, swap their positions. This picture, a typical example of which is 
displayed in Fig. 4, is quite similar to the classical description of collisions 
between solitons in the Korteweg - de Vries equation |22| . 

Collisions between pulses with opposite polarities and close velocities are 
found to be more inelastic in comparison with the unipolar pulses. In this 
case, dispersive wave packets with a considerable amplitude are generated, 
see an example in Fig. 5. Moreover, a new small-amplitude pulse with the 
negative polarity is also generated by the inelastic collision shown in Fig. 5. 

Another noteworthy feature revealed by the simulations is that the results 
of the collisions are very different depending on which solitons (faster or 
slower ones) have initially positive and negative polarities. This feature is 
obvious from the comparison of Figs. 5 and 6, which differ by the polarity 
reversal of the initial state. As is seen, in the latter case the inelasticity is 
much weaker, and, in particular, no additional pulse is generated. 

To further illustrate the inelasticity of the collisions, in Fig. 7 we display 
the shapes of the fields before and after collisions that were shown in Figs. 
5 and 6. In the former case (Fig. 7a), the amplitudes of the positive- and 
negative-polarity pulses decreases and increases, respectively, as a result of 
the collision. In the latter case (Fig. 7b), the result is the opposite. From 
here, we conclude that, in the course of the collision, the energy is always 
transferred to a more powerful pulse (the one with a larger value of a) . 

We also explored interactions between two pulses with equal values of a 
(hence, they have zero relative velocity), initially placed at various distances 
from each other. In Fig. 8, solid lines show trajectories of the motion of centers 
of the interacting pulses with different polarities, and dotted lines show the 
same for interacting unipolar pulses. In the case shown in this figure, the 
initial distance between the pulses is three half-widths of a steady-state pulse. 
At such a distance, tails of the pulses overlap considerably, giving rise to 
energy transfer from one pulse to the other (irrespective of the polarities, 
no visible interaction between pulses takes place if the initial distances is 
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increased to five half-widths). As a result, the pulses cease to be identical 
and start to separate (the pulse with a larger amplitude is moving slower). 
Eventually, the pulses assume permanent (but different) shapes after the 
separation. 

In the case of two unipolar pulses, energy is gained by the rear one, there- 
fore the trajectories of the pulses do not intersect in this case, see Fig. 8. 
On the contrary, in the case of opposite polarities the front positive-polarity 
pulse gains energy as a result of the interaction, while the absolute value of 
the amplitude of the negative-polarity rear pulse decreases and its velocity 
increases, hence the trajectories of the two pulses intersect. After their col- 
lision, there appears a new small-amplitude pulse with negative polarity (its 
trajectory is not shown in Fig. 8), cf. Figs. 5 and 7a. 

We stress that we have never observed formation of a bound state of 
pulses with equal or opposite polarities as a result of their interaction, nor 
was it possible to find bound states in any other way. 

To conclude this section, it is relevant to compare the present model with 
some others. Indeed, after transforming the system into the single equation 
(H), replacing ( and r by t and x, q by — u, and setting /i = one obtains 

dt dx dt dtdx 2 



which resembles to the so-called regularized long-wave equation |24], 



du du ^ du d 3 u 

dt dx dx dtdx 2 
While Eqs. ( |5B| ) and ([£|) have similar traveling-pulse solutions, the equations 
are not equivalent, hence they may produce very different dynamical effects. 



6 CONCLUSION 

We have introduced and analyzed a model for the propagation of ultrashort 
unipolar pulses of electromagnetic field in a material medium represented by 
anharmonic oscillators with quadratic and cubic nonlinearities. Two fami- 
lies of exact analytical solutions, with positive and negative polarities, have 
been found for moving solitary pulses. Direct simulations have demonstrated 
strong stability of the pulses against various perturbations. Collisions be- 
tween the pulses were also simulated in detail, showing that they interact 
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nearly elastically, irrespective of the pulses' relative polarity, unless their rel- 
ative velocity is very small. If the case of the small relative velocity, collisions 
are inelastic, resulting in generation of radiation and a new small-amplitude 
soliton. 
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7 FIGURE CAPTIONS 



Fig. 1. Numerically simulated propagation of the positive-polarity (a) and 
negative-polarity (b) ultrashort pulses predicted by the analytical solutions 
(H) and (|8D, © with a = 7.25 and fi = 2/9. 
Fig. 2. Numerically simulated evolution of the pulse with a = 5 and 
fi = 1/3, with a superimposed perturbation in the form of a quasi-harmonic 
wave packet. 

Fig. 3. The same as in Fig. 2, with an initial perturbation in the form of 
an additional rectangular pulse. 

Fig. 4. An example of a numerically simulated collision between two 
unipolar pulses with close velocities, corresponding to a = 5 and a = 4.61, 
at fi = 1/3. 

Fig. 5. Collision between pulses with the negative and positive polarities 
in the case of /i = 1/3. The pulse of negative polarity relates to a = 5 and 
pulse with positive polarity corresponds with a = 3.89 . 

Fig. 6. Collision of pulses with the same values of a and /i as in a fig. 5, 
but with reversed polarities. 

Fig. 7. The shapes of the pulses before and after the collisions, the panels 
(a) and (b) pertaining to the two cases shown in Figs. 5 and 6. 

Fig. 8. Trajectories of the motion of two interacting soliton with exactly 
equal initial velocities (equal values of a). The solid and dotted pairs of lines 
pertain, respectively, to the solitons with opposite and equal polarities. 
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